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Abstract 

We have calculated the perturbative corrections to all the structure functions in the 
semileptonic decays of a heavy quark. Assuming an arbitrary gluon mass as a technical 
tool allowed to obtain in parallel all the BLM corrections. We report the basic applica- 
tions, viz. perturbative corrections to the hadronic mass and energy moments with full 
dependence on the charge lepton energy cut. In the adopted scheme with the OPE mo- 
mentum scale separation around 1 GeV the perturbative corrections to {M%) are small 
and practically independent of E cnt ; the BLM corrections are small, too. The correc- 
tions to the second mass squared moment show some decrease with E cut consistent with 
the effect of the Darwin operator, within the previously estimated theoretical uncertainty. 
Perturbative corrections in the pole-type schemes appear significant and vary with Eg, de- 
creasing the moments at higher cuts. The hardness of hadronic moments is quantitatively 
illustrated for different cuts on E £ . 
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1 Introduction 



Inclusive decays of the heavy flavor hadrons are one of the most developed dynamic uses 
of the fundamental QCD in the short-distance regime. Application of the OPE allowed to 
address them at the nonperturbative level through the expansion in the inverse mass of 
the decaying quark: absence of the potentially largest 0( QCP ) corrections to the decay 
rates was established [1] and the leading power corrections to inclusive decay distributions 
were obtained through the expectation values of the local heavy quark operators in the 
decaying hadron [2-4]. As it has been known from the early days of QCD [5], the nonper- 
turbative treatment is best applied where usual perturbative corrections are comparable, 
or subdominant to the leading nonperturbative effects. This environment is often realized 
in the beauty decays once the OPE in the Wilsonian implementation is applied, isolating 
effects of large distances even from the traditional 'perturbative' corrections. 

High precision checks of theory and practical use of the inclusive decay distributions 
are now performed with new data of better reliability and of qualitatively different statis- 
tics. Matching this progress requires refinement of the theoretical predictions. Pertur- 
bative corrections to a number of decay distributions, first of all charge lepton spectrum 
in semileptonic decays of b quarks, have been calculated long ago [6], along with many 
other characteristics [7-9]. Most of the perturbative calculations, although often rather 
sophisticated, however aimed at a particular observable (total width, semileptonic spec- 
trum or its moments, angular asymmetries, etc.). As a result, in spite of the significant 
intellectual effort invested, they often remained highly specialized, without the possibility 
to be applied to only slightly modified observables. 

The hadronic mass and energy distributions of the final state in the semileptonic decays 
appeared very promising for scrutinizing nonperturbative QCD in the heavy quark system. 
In practice, experiments typically have to apply lower cuts on the energy of the charged 
lepton, to suppress backgrounds; this complicates the kinematics. Since, following the 
first papers on the dynamic heavy quark expansion, the power corrections were obtained 
directly for the semileptonic decay structure functions [4,10], calculating nonperturbative 
effects for any inclusive moment became straightforward, more or less regardless of the 
lepton energy cut. 

The perturbative effects, on the other hand, appeared in the role of the weak link - the 
perturbative corrections fully incorporating the cut in the lepton energy have generally 
been unknown for hadronic moments. 

To eliminate such an obstacle once and forever, we have calculated the perturbative 
corrections directly to the semileptonic decay structure functions. This allows straight- 
forward evaluation of the perturbative correction to all possible inclusive semileptonic 
distributions. Moreover, we allowed for an arbitrary fictitious gluon mass in the calcula- 
tions; this opens the possibility to calculate BLM corrections to an arbitrary order, or even 
to perform complete BLM resummation similar to the one accomplished in Refs. [11, 12], 
on the parallel footing. To be most universal for possible applications, following Ref. [4] 
we compute separately the structure functions induced by the vector and the axial- vector 
weak current, as well as their interference (w 3 ). Likewise all five structure functions are 
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calculated, so the results can be readily applied to the semileptonic decays into r-leptons, 
etc. The explicit structure functions can also be used for computing possible averages 
with the weight varying depending on the lepton energy cut. This is one of the possibili- 
ties to optimize the trade-off between the sensitivity to the heavy quark parameters and 
suppressing experimental backgrounds. 

When this study was in the completing phase, the preprint by M. Trott [13] appeared, 
where the first results on the perturbative structure functions proper were presented. This 
is clearly an important step in completing the theoretical toolbox required for precision 
analysis of B decays. While addressing in general closely related subjects, we differ with 
Ref. [13] in a number of points calculation- wise, and also in our applications. Aiming for 
ultimate flexibility in future applications, we have calculated all five structure functions 
separately for the vector and the axial-vector currents. Moreover, our results are equally 
applicable to the first-order perturbative corrections and to all BLM corrections. 

Although having independent analytic calculations is important for cross-checks, at 
this point we have not attempted to compare our results with those reported in Ref. [13], 
even for the pure one- loop limit. Programming alternative expressions anew for numerical 
evaluation would represents a time-consuming process; it is not easily safeguarded against 
possible typos at this technical, yet unavoidable step. To check our calculations we had 
used other opportunities mentioned in Sect. 5, made available through our previous studies 
of perturbative and power corrections in the semileptonic decays of heavy flavors, or 
provided by the OPE applied to perturbation theory. 

This paper reports a certain progress in the project presently carried out together with 
Paolo Gambino, dedicated to the precision analysis of the inclusive B decays. It uses the 
Wilsonian implementation of the OPE to combine perturbative and power-suppressed 
effects. In this paper we mostly address the qualitative features and report only the basic 
applications important for the ongoing experimental analyses, which complement our 
recent publication [14]. 1 The comprehensive presentation is planned for the forthcoming 
paper [16]. 

2 Inclusive B decays in the perturbative expansion 

The semileptonic structure functions are defined as the absorptive part of the covariant 
structures appearing in the decomposition of the forward scattering amplitude of the two 
weak currents off the B meson: 

M<7 2 ,<7o) = ^<£| Jd 4 xe-^iT{j,(x),Jt(0)}\B), (1) 

where J M is generally the cry^fr or/and cq^b current (b — >u decay simply corresponds to 
m c ^0). Following the standard notations of Ref. [4] we put 

/V = -hi9/w + h 2 v fM v u - ih 3 e^ a(3 v a q p + h^q v + h^Vy + v^q u ) , (2) 

1 Om: notations follow that paper; for the nonperturbative operators we consistently use notations of 
Ref. [15]. 
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and 

Wi(q 2 , g ) = 2 lmhi(q 2 , g ) , (3) 

with denoting the 4- velocity of the decaying meson (or quark, for perturbative calcula- 
tions); q 2 and go have the meaning of the invariant mass squared and combined energy of 
the lepton pair, respectively All the decay distributions with light leptons are expressed, 
for example in terms of the first three structure functions [4]: 



d 3 r _ G 2 F \V cb \ 



dE e dq 2 dq 32tt 4 



#(q -E e -£r)#(E e )#(q 2 ) x 

{2gV + [4E e (q -E e )-q 2 ]w 2 + 2q 2 (2E £ -q )w 3 } ; (4) 



the total integrated width without cuts on lepton energy depends only on w\ and W2'- 
„ G%\Vrh\ 2 r m l , , f m "> , rr, — -/, , „ , 1 



'F\ v cb\ 

16tt 4 



rmt rnif, i / 1 \ 

/ dg 2 / ^-dgoygo-g 2 g 2 wi(g 2 ,g ) + -(go-g 2 )w 2 (g 2 ,go) • (5) 

JO Jqo>Jq 2 V 3 / 



Perturbative structure functions to order a\ consist of virtual corrections oc 5(qo — 

b 2mb — -), and bremsstrahlung contributions with g < — — > where A is the 

gluon mass. In the limit A 2 — > virtual corrections logarithmically diverge. The divergent 
part is given precisely by the tree-level (free quark) structure functions, see Ref. [4], 
Eqs. (Al-All), with the universal coefficient depending on g 2 . The real emission part at 
small A 2 has a singularity near the free quark kinematics, 

^(g 2 ,g )oc _ J_ +g2 _ mg , (6) 



go 



2m b 



with the commensurate coefficient, so that the divergence cancels once integration over 
g is performed. At g 2 approaching {m h —m c ) 2 the coefficient vanishes, therefore shrinking 
domain of integration over g does not affect the cancellation of divergences. 

Perturbative structure functions themselves are strongly infrared-sensitive in the dom- 
inant domain close to the free quark decay kinematics. However, the integrals entering the 
experimentally measured inclusive moments are not. Moreover, their infrared sensitivity 
is governed by the OPE, and it is greatly reduced when using the Wilsonian separation 
between perturbative and power corrections based on the momentum scale a particular 
contribution originates from. This separation is in practice done directly for the moments 
rather than for the structure functions themselves. 

The analytic expressions for the 0(a s ) bremsstrahlung structure functions Wk(q 2 ,qo) 
are relatively simple containing at worst logs of the kinematics-related square roots, even 
at arbitrary gluon mass. Yet they are lengthy for separate structure functions, consisting 
of many terms at A 2 7^ 0. They simplify significantly at A 2 = 0. Virtual corrections at 
A 2 — > are well known one loop renormalization of quark currents containing dilogs at 
worst. At arbitrary A 2 , however the one- loop vertices with three different internal masses 
and a general momentum transfer are too special functions. In practice, we represent 
them as one-dimensional integrals over a single Feynman parameter u, with 0<u<l. 
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Obtaining a general moment with the lower cut on lepton energy we, therefore integrate 
the explicit structure functions over q and q 2 (for bremsstrahlung) , or over u and q 2 (for 
virtual corrections), with the weight which generally is a polynomial in (m b — go) an d 
(m 2 + q 2 — 2m b qo — m 2 ). The lepton cut enters through the concrete weight following from 
Eq. (4): it is obtained by integrating the coefficients from E[ to E e + , where 

Ei = q J^MEl , E[ = max{ £ cut , E _ } . (7) 
the weights are simple polynomials in E[ or E l + . 

3 Applications 

Recent experiments at B factories provide data of impressive precision and quality, which 
allows for stringent tests of the heavy quark expansion for inclusive decays, direct exper- 
imental extraction of many heavy quark parameters and for robust defendable extraction 
of V c b and V ub from the integrated decays rates. The current precision requires to fully 
implement the applied cut on the charge lepton energy in the theoretical calculations. 
In the recent publication [14] the OPE-based predictions were given for various moments 
with cuts, in the framework which does not rely on assuming charm to be a heavy quark, 
but only expanding in l/m b . Perturbation theory- wise we used the Wilsonian approach 
which assumes excluding soft gluon contributions from the coefficient functions. This 
rendered the perturbative corrections small in size and presumably stable against higher- 
order corrections. 

For hadronic mass moments Ref. [14] evaluated the perturbative corrections without 
the cut on lepton energy, since the full perturbative corrections with cuts had not been 
known. Although some simpler parts of the corrections had been calculated, including 
them would not be consistent: in hadronic moments the terms proportional to different 
powers of Mg-mj essentially mix under renormalization. In particular, for the average 
hadronic mass square, 

(M 2 X ) =m 2 c + (M B -m b ) 2 + 2(M B -m b )(E x ) + (m 2 x -m 2 c ) (8) 

the last term having the meaning of the excess of the combined parton invariant mass at 
the quark level over ml, is fed under renormalization by the preceding term driven by 
(E x ), the average hadron energy at the parton level [14]. 

The justification for approximation neglecting E e -cut in the perturbative corrections 
was, again, using the Wilsonian version of the OPE. The absolute size of the perturbative 
corrections is suppressed here, so even their noticeable variation with the cut was not 
expected to produce a significant bias. Moreover, it is the soft parton processes that are 
most sensitive to the particular kinematics. The truly hard gluon effects contributing the 
perturbative corrections in our approach, are expected to be less dependent on the details 
of the kinematics as long as the process in question remains sufficiently 'hard'. 
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Figure 1: Dependence of the perturbative corrections to (M|-) on E cut , for \i = 1 GeV (red 
curve). Long-dashed and short-dashed curves show the separate contributions of A c and E x , 
respectively. Blue line (top) gives the perturbative correction in the pole-type scheme. 

With full perturbative expressions available, we can check validity of those tempo- 
rary assumptions. We indeed found that in our scheme the perturbative corrections to 
(Mf) are practically independent of E cnt in the whole domain E cnt ^ 1.4 GeV. A more 
noticeable - yet still not very significant - variation with E cut is observed for the second 
hadronic moment with respect to average, ((Mf- — (M^)) 2 ). This is expected, since this 
higher moment is more sensitive to the Darwin operator, for which we did not actually 
remove the corresponding soft contributions from the leading-order coefficient function. 
The observed -E cut -dependence, in particular the sign and the size are compatible with 
such a contribution from the Darwin expectation value: its coefficient function is negative 
and increases in magnitude for larger E cut [14], see, e.g. Table 6 of Ref. [14]. 

The full corrections are illustrated in Figs. 1 and 2 as functions of E cnt . For comparison 
we also show the corresponding effect in the 'usual' (pole-type) perturbative corrections 
without Wilsonian separation. It is evident that for the first moment the effect is dra- 
matically different. 

The third moment of the hadronic invariant mass squared places, for the perturbative 
corrections, most weight on harder gluons. In its soft part it is mostly sensitive to the 
D = 6 Darwin operator. Therefore, using the Wilsonian prescription only for the effects 
scaling like l/m,b and 1/m 2 does not make a noticeable difference, Fig. 3. The formal 
perturbative contribution is then significantly dependent on the cut; 2 one should keep 
in mind, though that a significant fraction of it still comes from gluon momenta below 
1 GeV. (For illustration we show the result of subtracting the contribution of the Darwin 
expectation value of 0.062 GeV 3 corresponding to /i = 0.9 GeV in fixed-order perturbation 
theory). In practice, the third moment so far is important mainly as an estimate of the 

2 Ref. [14] evaluating the third hadronic mass moment did not include the perturbative corrections 
to the terms oc (Mb — mf,) fc with k > 1, even at zero cut. We have computed them, and numerically 
they turned out to be about 1 GeV 6 (the dashed lines excluding the highest one), decreasing fast from 
E cut > 500 MeV. 
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Figure 2: Perturbative corrections to ([JW| — (M|-)] 2 ), for /x = 1 GeV (red curve), and their 
breakdown showing the separate contributions oc (Mb— mi)) k , k = 0, 1, 2 (dashed curves). Blue 
line (top) refers to the pole scheme. Orange dashed-dotted line shows subtracting the soft piece 
of 0.062 GeV 3 of the Darwin expectation value; it would correspond to Wilsonian pp. normalized 



scale of p 3 D ; as had been pointed out [17,18] other unaccounted effects introduce theoretical 
uncertainties of the same magnitude. 

4 Running ct s and BLM corrections 

We are also in the position to calculate the so-called BLM corrections - the effects ac- 
counting for running of the strong coupling a s in one-loop perturbative diagrams. Such 
effects are typically quite significant where the perturbative calculations do not eliminate 
explicitly the low-momentum domain; however, they are moderate, or even small in size 
in the appropriate Wilsonian scheme. A dedicated discussion for the case of the total 
semileptonic width was given in Refs. [11,12]. We leave a detailed implementation of this 
method for the future publication [16], and here place the emphasis on the qualitative 
features. They allow to assess the possible accuracy of the theoretical predictions without 
much numerology. 

The technique for the BLM summation in the Wilsonian approach was first discussed 
in Ref. [19]; it is concisely presented in recent Refs. [11,12], and we do not describe it 
here. Neither we detail numerical predictions or address all experimentally interesting 
observables. We also largely leave aside more subtle theoretical aspects related to full 
BLM resummation for higher moments with only a few lowest nonperturbative operators 
included. Our main interest lies in the most sensitive case for theory, viz. the width and 
the first hadronic moment with the lepton energy cut. 

We start with the decay width itself when the lepton energy cut is imposed. A con- 
venient way to visualize the contribution of various momentum scales Q 2 in one-loop 



at 0.9 GeV. 
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Figure 3: Similar plot for the third invariant mass moment ([M^ — (M^)] 3 }. Total perturbative 
corrections are reduced (orange dashed-dotted line) using p|>(0.9GeV) compared to the case of 
fP D (red solid line). Blue solid line corresponds to the pole scheme. 

perturbation theory for an observable A is provide by the distribution 

dlnQ 2 ~ din A 2 

where A^^X 2 ) is the first-order perturbative coefficient calculated with the fictitious non- 
zero gluon mass [20,21]. This representation has some limitations, yet is quite suitable 
for qualitative purposes and in adopted in the plots presented below. 

Fig. 4 shows the distribution over the gluon virtualities for the two extreme cases, 
E C ut = (no cut) and E cut = 1.5 GeV. Without the scale separation, dashed lines, the 
perturbative corrections are clearly not too well behaved, since the major contribution 
comes from the gluon momenta noticeably below 1 GeV. This is an artefact brought in 
by the pole masses [22,23] used by conventional perturbative diagrams. Fig. 4 shows that 
applying Wilsonian procedure quite effectively eliminates this domain (solid line). 

At first glance, there is no much difference between the perturbative corrections for 
-Ecut = 1-5 GeV compared to the total width without kinematic restrictions. However, 
the low-momentum tail is much higher for T s \(E e > 1.5 GeV). This leads to a far more 
significant impact of the soft physics, the fact deduced independently from the growth of 
the Wilson coefficient for the higher-dimension nonperturbative operators (e.g., Darwin). 
As anticipated, the effect is more significant for moments of the distributions. We will 
now look closer at the average hadron invariant mass squared (M|). 

In the OPE, M\ consists of two distinct dynamic pieces (they mix under the renor- 
malization of the effective low-scale QCD for heavy quark): 

(M 2 X ) = m 2 c + (M B -m b ) 2 + {M B -m b )(2E x ) + (m 2 x -m 2 c ) (10) 
E x = m b -q , A c = m 2 x -m\ = m 2 b + q 2 -2m b q -m 2 c . 



X 2 -- 



(9) 
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Figure 4: Distribution over the gluon momentum for the total width without a cut (blue) and 
for E CVLt = 1.5 GeV (red). Solid lines correspond to the cutoff scale (jl = 1 GeV, dashed lines show 
the unsuppressed contributions in the pole-type schemes. The area bounded by a curve yields 
the overall first-order perturbative coefficient. 

The first one, (MB — mb)(2E x ) dominates corrections to (M|) [24]. However, among the 
perturbative corrections those to (A c ) appear most significant. Therefore, we start with 
the second piece, (A c ). 

At zero lepton cut, A c which in the conventional perturbative diagrams is described by 
real gluon emission, is very soft, see Fig. 5: the four-body phase space for b^cg+iv is very 
sensitive to the gluon mass even when the latter is only a few hundred MeV. Presence of 
such corrections would be a disaster for precision calculations. However, in QCD the phase 
space-unsuppressed soft gluon emissions feeding m 2 x are related by gauge symmetry to the 
Coulomb self-energy of the heavy quark [15]. The Wilsonian treatment then effectively 
eliminates the infrared domain, as illustrated by the solid curve corresponding to the 
separation scale \x = 0.9 GeV. 

The similar distribution for the cut at 1.5 GeV is shown in Fig. 6. Applying the same 
Wilsonian cutoff at 0.9 GeV we eliminate the deep infrared domain as well, yet paying the 
price of quite unphysical spike in the contribution of gluons at the scale around 1 GeV. 
It disappears only when \x is lowered down to 0.55 GeV. This behavior does not mean 
that the masses and higher- dimension operators normalized at, say 1 GeV cannot be used 
for perturbative calculations. The perturbative result can be expressed in terms of any 
short-distance masses, for instance m\,{m^) and m c (m&) corresponding to /i in a few GeV 
range. It rather demonstrates that the effect of the domain of momenta between 600 MeV 
and 1 GeV is not fairly described by only the leading 1/mb and 1 / ml contributions - the 
terms scaling like higher powers of l/m^ are comparable, or even dominate them. This 
is a direct consequence of the deteriorating hardness of the inclusive width with the high 
cut, which effectively introduces another, much lower than mj mass scale parameter in 
the OPE [17,18,25]. 

For milder cuts with lower cutoff energy the convergence of the power expansion 
improves; for instance, Fig. 7 shows the similar distribution at E cnt = 0.9 GeV where it 
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Figure 5: Similar distribution for the domi- 
nant perturbative contribution A c to (M^), 
without lepton energy cut. Solid line is for 
/i = 0.9 GeV, dashed for fi = (pole scheme). 




Figure 6: The case of E cnt = 1.5 GeV. Dot- 
ted line refers to the same separation scale 
// = 0.9 GeV; solid line shows the choice of 
^ = 0.55 GeV. 




-1 . 5 



0.5 



mi. 



Figure 7: The case of a mild cut at £^ = 0.9 GeV, using fj,= 0.8 GeV. 



still looks perfect at /x = 0.8 GeV. 

The similar scrutiny can be applied to another short-distance component of (M|), the 
perturbative corrections to (2E X ). Its naive 0(a s ) evaluation turns out quite small and 
might be thought to be of no practical interest. However, we see from Figs. 8 and 9 that 
this is a result of cancellations between softer and harder gluons; as such it is vitiated 
already when running of a s is accounted for. The effect of the higher-order corrections, 
taken at face value appears dramatic here. However, in the Wilsonian approach all these 
corrections can be readily accounted for, and lead only to minute changes in (M|) for 
reasonably placed cuts on E e . 

Similar in spirit and technically anatomy can, in principle, be applied to higher 
hadronic moments with cuts as well, where the terms accompanying different powers 
of (M^— ra&(//)) should be analyzed separately (say, for the second mass squared moment 
there are three terms with the power from to 2). We have to bear in mind, however 
that the higher moments more seriously depend on the Darwin operator. Therefore, ac- 
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Figure 8: The gluon virtuality distribution Figure 9: The same for E cut = 1.5 GeV, ap- 

for (E x ) without lepton energy cut; separa- plying ^ = 0.45 GeV. 

tion scale ^ = 0.7 GeV. 



counting for running of a s would not be too meaningful without extending the Wilsonian 
treatment to higher- dimension operators, first of all to the Darwin operator. We will 
present the corresponding results in the forthcoming publication [16]. 
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Figure 10: Effect of the BLM corrections on (Mj^). Upper curves show the combined contri- 
bution: solid for the fixed-order a s = 0.3 evaluation, dashed line for order O((3oa 2 ) and dotted 
line at 0(f3 2 o? s ). Lower curves are similar effects for the contribution oc (Mb — mi,) alone. 

Here we illustrate the numerical results for (M|). Fig. 10 shows the breakdown of 
the predictions for the perturbative corrections to (M|) as a function of the cut on lep- 
ton energy, using the fixed-order perturbative estimates with a s = 0.3 [14], a generally 
appropriate choice for beauty decays. For comparison we show also the predictions ob- 
tained by combining the first-order 0(a s ) term evaluated with admittedly too low a value 
a^ s = 0.22, with the second-order BLM correction, O((3oa 2 ): 



(M 2 X ) = m 2 c (fi) + (M B -m b (fi)) 2 

+ (M B -m b (fi))m b lE{ 0) + C F 



+ m 2 h C F 



^ + f(-)^ 2) + (f) 2 (- 



3 £f 3) 



+ 



.(3) 
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a s = af s (m b ) , 



C F 



4 
3' 



A) 
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Such an approximation is routinely applied in lieu of a more thoughtful choice of a s 
where the corresponding BLM correction is known. We can even add the higher-order 
BLM terms (dotted lines in Fig. 10). 

It is clear that the higher-order corrections are under good control here. Taking the 
results literally, the effective value of a s to be used in the fixed-order calculations is only 
slightly larger than 0.3 slowly increasing at higher cuts, and does not exceed 0.38, well 
within the interval allowed for in Ref. [14]. 

To have an unbiased comparison, however one should recall that the part of the shift 
when using the BLM-improved form is actually offset by the change in the commensurate 
value of the effective Darwin expectation value: 



Pi 



' p^(lGeV) - 0.085 GeV 3 

p^(l GeV) -0.12 GeV 3 
k p%{\ GeV) -0.16 GeV 3 



fixed order, 
O((3 a 2 s ), 



o s = 0.3 
a, = 0.22 
a, = 0.22 



(12) 



For the first-order BLM improvement this shift lies below theoretical accuracy and we 
usually neglect it, but it is relevant when assessing the net effect of the BLM improvement. 
For higher orders it becomes significant. 



Table I. Perturbative coefficients for (M|-) 



at m c = 1.16GeV, m b = 4.6GeV, /u = lGeV 



Scut, GeV 


A? 


Ai 2) 


A? 


Ai 4) 




E? 


E™ 




E^ 





0.021 


0.036 


0.065 


0.3 


0.846 


0.0098 


0.061 


0.097 


0.49 


0.6 


0.02 


0.035 


0.064 


0.3 


0.839 


0.0085 


0.058 


0.095 


0.49 


0.9 


0.019 


0.034 


0.063 


0.3 


0.827 


0.0069 


0.056 


0.095 


0.51 


1.2 


0.018 


0.036 


0.067 


0.33 


0.811 


0.0069 


0.06 


0.11 


0.58 


1.5 


0.023 


0.046 


0.087 


0.44 


0.796 


0.015 


0.088 


0.15 


0.8 



Relegating the detailed analysis of other inclusive semileptonic averages including 
higher hadronic moments, to the dedicated publication [16], here we give a few perturba- 
tive coefficients, including BLM corrections, for the first, Table 1, and second, Tables 2a 
and 2b, hadronic mass square moments, at different lepton energy cut. The coefficients 
for the second moment are defined in analogy to the case of the first hadronic moment, 
Eq. (12), however they proliferate: 



+ 



{[M 2 X - (Ml)] 2 ) = (M B -m b ) 2 m 2 b \E^ +C F 



•-E. 



(i) 



§0 
2 



+ (M B -m b )m b C F 



lA {1) 
^11 



+ 



(-) Z Ai? + (fH t )"A?> 



+ 



m b C F 



+ f fr) 2 Af + (f) 2 (-) 3 Af + ... 



(13) 
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It should be stressed that higher-order coefficients are quoted only for illustration, and 
should not really be added to the corrections without carefully adjusting for the effect of 
the residual infrared piece. 

The estimated perturbative corrections to the second ((M|- (M|)) 2 ) and to the third 
((M|— (M|)) 3 ) hadronic mass moments depending on the lepton cut, and their breakdown 
into separate terms are shown in Figs. 2 and 3, both in fixed-order perturbation theory. 
We also show the results in the pole scheme (/x = 0) assuming, however, the same values 
of the quark masses. For numerical evaluations to order a] we commonly use a s = 0.3, 
and assume the values of the short- distance heavy quark masses 

m c (lGeV) = 1.16 GeV, m b {\ GeV) = 4.60 GeV 

preferred by experiment. 



Table 2a. Perturbative coefficients for ([M| - (M^)] 2 ), in the same setting 



-^cut 


A^ 


A?> 


A?> 




A (D 


A(2) 
^11 


A(3) 


A(4) 
^11 





-0.00082 


-0.0035 


-0.0077 


-0.046 


-0.0023 


-0.01 


-0.023 


-0.14 


0.6 


-0.0016 


-0.0044 


-0.0091 


-0.052 


-0.0043 


-0.013 


-0.027 


-0.16 


0.9 


-0.0026 


-0.0057 


-0.011 


-0.06 


-0.007 


-0.016 


-0.033 


-0.18 


1.2 


-0.0036 


-0.0073 


-0.014 


-0.073 


-0.0099 


-0.02 


-0.039 


-0.21 


1.5 


-0.0046 


-0.0091 


-0.017 


-0.089 


-0.012 


-0.024 


-0.046 


-0.24 



Table 2b. The coefficients for the term oc {Mb— mi)) 2 



E C ut 


4 0) 


E? 


4 2) 


4 3) 


4 4) 





0.194 


-0.0038 


-0.017 


-0.029 


-0.17 


0.6 


0.2 


-0.0054 


-0.02 


-0.032 


-0.17 


0.9 


0.208 


-0.0076 


-0.024 


-0.035 


-0.18 


1.2 


0.218 


-0.01 


-0.03 


-0.041 


-0.2 


1.5 


0.226 


-0.013 


-0.035 


-0.047 


-0.23 



The first-order perturbative corrections to the second mass squared moment happen 
to be rather suppressed near zero cut on E e . To some extent this is accidental and takes 
place just for our choice of /x around 1 GeV: it would be nearly complete for a somewhat 
smaller /j, and softer for a larger separation scale. The BLM corrections at small E cut 
then seem to have a large relative impact, nearly doubling the perturbative correction. 
However, as seen from Fig. 2, the cancellation fades out at a relatively low cut, and the 
significance of the BLM corrections becomes moderate already at E cut > 600 MeV. Above 
0.9 GeV the face value of the effective a s for the one-loop contribution is about 0.35 - even 
discarding the possible effect of redefining the appropriate Darwin expectation value p 3 D . 
Therefore, it seems likely that the actual perturbative corrections to the second moment 
is even flatter than follows from Fig. 2. 
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It has been pointed out in Ref. [14] that, in view of possible cancellations in the first- 
order perturbative coefficients, it is unsafe to estimate the uncertainty associated with 
uncalculated perturbative corrections simply as a fixed fraction of the one- loop result. 
This may particularly apply to the schemes other than the pole one. Ref. [14] suggested 
that an additional uncertainty should be added obtained examining certain 'perturbative' 
variations of the nonperturbative expectation values, to evade such special cases. We 
see that this recipe works well here: while even a 50% variation in the perturbative 
corrections to the second mass moment near zero cut might easily underestimate their 
actual uncertainty, considering the effect of the Darwin operator gives the right estimate 
of the size, in fact on the safer side. 

5 Discussions 

We have calculated the one-loop perturbative corrections to all the semileptonic decay 
structure functions of the heavy quarks, in the form allowing to obtain all BLM correc- 
tions on the same footing. This parallels the work done in Ref. [26] almost a decade 
ago which, however was applicable only to the total semileptonic width. With the com- 
plete expressions for the power corrections to the structure functions through order 1 / TOq 
available [4,10], there remain no practical obstacles in calculating all sufficiently inclusive 
semileptonic distributions with significant precision. As has been emphasized in Ref. [12] 
in respect to precision determination of \V c b\, now the real limitation in many cases be- 
comes perturbative corrections to the Wilson coefficients of the power-suppressed nonper- 
turbative operators; sometimes the effects of the nonperturbative four-fermion averages 
with charm fields (called 'intrinsic charm' in Ref. [12]) may be noticeable. 

We think that the program of precision extraction of higher nonperturbative heavy 
quark parameters of B mesons will be on agenda regardless of calculation the perturbative 
corrections with all kinematic constraints - in many instances already power corrections 
proper lead to significant uncertainties. In this respect studying at S-factories the mod- 
ified higher hadronic mass-energy moments M x detailed in Ref. [14] looks promising. 
From the perturbative point of view they are similar, or even simpler than the standard 
moments of M x , and our analysis applies to them in full. 

We plan to present both perturbative and nonperturbative corrections in an (possibly 
cumbersome) analytic form in the forthcoming publication [16], and to simultaneously 
provide ready-to-use computer programs making numerical evaluations an automated 
procedure. 

Calculation of the perturbative 6-decay structure functions to order a] is rather 
straightforward if only lengthy. The main problem is to control possible typos, in particu- 
lar in the computer program which would evaluate the corrections. A single omission can 
both radically change the emerging numerical result, or to be nearly invisible numerically 
in a particular kinematic setting. To safeguard our results against such complications, we 
have applied a number of checks. 

Firstly, a regular limit for a massless gluon was verified for all the moments, and the 
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logarithmic dependence of separately the virtual and the bremsstrahlung contributions at 
given q 2 was checked to match the known classical bremsstrahlung radiation divergence. 
Furthermore, the limit A 2 — > in the sum of the contributions for a particular moment, 
is slow containing the terms oca/A 2 , rather than only those scaling like A 2 . These terms, 
however are controlled by the OPE - they reside solely in the corresponding dependence of 
the pole quark masses on A 2 . Moreover, the OPE ensures that in the first-loop corrections 
terms oc A 2 In A 2 do not appear. This provided the serious test for the moments - once 
passing to the Wilsonian scheme at a given gluon mass (and such a translation involves 
coefficients with a nontrivial kinematic dependence), we observed a straight linear in A 2 
behavior at the small gluon mass. 

At arbitrary gluon mass of order m c or below we observed that the total width we 
numerically calculate precisely coincides with the same width we used in Refs. [11,12]; 
those expressions were derived in Ref. [26] using, generally, the different set of integration 
variables. The total width, however depends only on w\ and w 2 , see Eq. (5). We also 
compared the numerical results of the perturbative corrections to the width with an 
arbitrary lepton energy cut at A 2 — > with the computations used previously [14] based 
on the well-established analytic expressions of Ref. [6] ; this verified W3 as well. Together 
with a number of more technical cross checks, a confidence has been gained that possible 
typos in the program for evaluating moments were all eliminated. 

Although perturbative effects are conceptually simple, getting accurate predictions 
requires thoughtful combining perturbative and nonperturbative corrections in the op- 
timal way. Our experience commencing with an early paper [27] suggests that using 
the literal Wilsonian prescription which separates various contributions based on their 
intrinsic momentum scale, is the way to get reliable predictions. In return, applying it 
often yields accurate results even with minimal computational efforts, or using simplifying 
approximations. 

In the present paper we have reported some results for the dependence of the per- 
turbative corrections to the hadronic mass moments, on the charge lepton energy cut; 
an extensive analysis is in preparation [16]. In accord with the expectations which con- 
sidered peculiarity of the adopted OPE implementation, we found weak variations, for 
the safe intervals in E cut , well within our estimates of the overall accuracy theory in its 
present form can realistically provide [14]. The perturbative corrections to (M|) actually 
turned out nearly constant, varying even much less than could be anticipated. The ^cut- 
variation in the second mass moment, ((M|- - (M^)) 2 ) is more noticeable, consistent in 
sign and magnitude with the cut-dependent contribution from the Darwin operator. We 
actually found indications that they may be further flatten when the BLM corrections are 
incorporated: the latter eliminate a somewhat accidental cancellation at very low E cut . 

It should be emphasized that such a moderate sensitivity resulted from eliminating soft 
gluon effects from the perturbative diagrams; in the schemes without such a separation 
(we generically refer to them as pole-type) the variation, as a rule, is far more pronounced. 
We also find that sensible BLM improvement of the one-loop estimates is possible within 
the Wilsonian approach, although it does not seem to have dramatic impact on the a 
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priori safe moments. Nevertheless, we point out that there are strong cancellations in 
the perturbative corrections to the hadron energy (E x ) between different domains of 
integration. They are vitiated once running of a s is accounted for, which strongly enhances 
their effect. Their significance have often been underestimated. 

Our predictions [14,25] both for the absolute values of the hadronic moments, and for 
their cut-dependence were in a qualitatively good agreement with the preliminary data 
reported by BaBar and CLEO. Having calculated the full perturbative corrections, we did 
not find effects which would be unexpectedly large or would show a surprising behavior; 
therefore, we expect the agreement to persist, or possibly even strengthen. More precise 
verifications of the theory including fits to many measured observables, is to follow from 
the dedicated experimental analyses. 
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